Conformal Behavior in High Temperature QCD 
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We give a new perspective on the properties of quarks and gluons at finite temperature T in 
Nf — 2 ~ 6 QCD. We point out the existence of an 1R fixed point for the gauge coupling constant 
at T > T c (T c is the chiral phase transition temperature). Based on this observation we predict 
theoretically and verify numerically that the correlation functions of a meson G.s(£) at T/T c > 1 
decay with a power-law corrected Yukawa-type decaying form, Gh(£) = c exp(— mHt)/t a in the 
"conformal region" defined by mu < cAir, where Air is the IR cutoff, ran is the ground state 
energy in the cannel H and c is a constant of order 1. The decaying form is the characteristics of 
conformal theories with an IR cutoff. We discuss in detail how the resulting hyperscaling relation of 
physical observables may modify the existing argument about the order of the chiral phase transition 
in the Nf = 2 case. 



The properties of quarks and gluons at high tempera- 
ture are key ingredients for understanding the evolution 
of the Universe and the heavy ion collision experiment. 
Lattice QCD is the most reliable formulation of QCD for 
the investigation of non-perturbative properties of quarks 
and gluons, and indeed a lot of works from the early stage 
of lattice gauge theories were devoted to clarify them [l[ . 
Although many interesting and useful results were ob- 
tained, still there remain unsolved issues. Among other 
things, it is a long standing issue to determine the or- 
der of the chiral transition for Nf = 2 QCD, which has 
recently attracted a lot of attention [2| . 

In this article we give a new perspective on the proper- 
ties of quark-gluon state above the chiral phase transition 
temperature T c , and thereby give a new insight into the 
order of the chiral phase transition. Our main idea is to 
use the concept of the conformal field theory with an IR 
cutoff. After presenting the renormalization group argu- 
ment and conjecturing that the propagator decays with 
a power-law corrected Yukawa-type decaying form, we 
discuss numerical results of lattice calculations and their 
physical implications on the chiral transition in Nf = 2 
QCD. Some preliminary results have been presented in 

We discuss lattice QCD at high temperature for small 
Nf (2 < Nf < 6) massless fermions in fundamental 
representation where the chiral phase transition occurs 
at some critical temperature T c . Our general argument 
that follows can be applied to any number of flavors 
(2 < Nf < 6) with any formulation of gauge theories 
on the lattice. 

For numerical simulations in this article we take Nf = 
2 and employ the Wilson quark action and the standard 
onc-plaquette gauge action on the Euclidean lattice of the 
size N x = N y — N z = N and N t , with the lattice spacing 
a. We impose an anti-periodic boundary condition in the 
time direction for fermion fields and periodic boundary 
conditions otherwise. In order to obtain physical quan- 
tities at temperature T, we have to take the thermody- 
namic limit N — > co and the continuum limit a — > 0, 



keeping Nt a = 1/T fixed. The order of the thermody- 
namic limit and the continuum limit can be exchanged. 

The theory is defined by two parameters; the bare cou- 
pling constant go and the bare degenerate quark mass too 
at ultraviolet (UV) cutoff. We also use, instead of go and 
to , (3 = 6/go and K = l/2(moa + 4). We define the 
quark mass m q through Ward-Takahashi identities 0, Q 
with renormalization constants being suppressed. 

We first consider the case that the renormalized quark 
mass is zero. Then the renormalized coupling constant is 
the only relevant variable in the theory. A running cou- 
pling constant g([i;T) at temperature T can be defined 
as in the case of T = 0. The following discussion can 
be applied to any definition of the running coupling con- 
stant g(fi;T). One promising way is the Wilson MCRG 
method in order to investigate the infrared behavior. 

In most of the existing literatures (see e.g. ref. @), the 
authors discussed the running coupling constant g(/j,;T) 
mainly in the UV regime. However, what we would like 
to understand is the IR behavior in connection with the 
issues raised in the introduction. 

In the UV regime, since the theory is asymptotically 
free, the running coupling constant at finite T can be 
expressed as a power series of the running coupling con- 
stant at T = as long as g is small [7| . The leading term 
is universal in the limit g — > 0. 

However, in the IR region, g(fi;T) is quite different 
from g(n\T = 0), since the IR cutoff Air in the time 
direction is T, while the IR cutoff is zero at zero temper- 
ature. Furthermore, when T/T c > 1, where the quark 
is not confined, the running coupling constant g(/x; T) 
cannot be arbitrarily large. This means that there is an 
IR fixed point with non-trivial zero of the beta function 
when T/T c > 1. This is the key observation in this arti- 
cle. 

As long as T < T c , the beta function is negative all 
through g. As the temperature is increased further, the 
form of the beta function will change as in FigQ] (left) 
When T > T c but T ~ T c , the beta function changes 
the sign from negative to positive at large g; (middle) 




gfaT) 




T»T, 



g(v,T) 




g(v;T) 



FIG. 1. The beta function f3(g(fj,;T)) at finite temperature T 



As temperature increases the fixed point moves toward 
smaller g; (right) When T ^> T c it changes the sign at 
small g. 

In Lorcntz invariant theories, the existence of an IR 
fixed point implies that the theory is conformal (see e.g. 
[8[ for a review). In our case with finite temperature, 
we argue that the consequence of the IR fixed point for 
T > T c naturally leads to the conformal field theories 
with an IR cutoff we have recently investigated in J10| . 
There, the dynamics of conformal theories was realized in 
the SU(JV) gauge theory, when the number of flavors Nf 
is within the conformal window Q , and the finite lattice 
size introduced the IR cutoff. The discussion presented 
there can be applied here with small modifications. We 
refer ref. |1(| f° r * ne KG argument which is a background 
of the predictions made in ref. [10( and here. 

To verify the above idea of conformal field theories with 
an IR cutoff realized in high temperature QCD, we will 
study the propagator of a local meson operator Gn(t) on 
a finite lattice for T > T c with T c defined for the chiral 
transition temperature in m q —> limit. The analogy 
to the conformal field theories with an IR cutoff defined 
for the large Nf (within the conformal window) QCD 
implies that the behavior of Gjj(t) qualitatively differs 
depending on whether the quark mass is smaller than 
the critical mass or not. 

When the theory is in the relatively heavy quark re- 
gion, it decays exponentially at large t as 



Now let us discuss the results of our numerical cal- 



G H it) = c H exp(-mfli), 



(1) 



where mu is the energy of the ground state in the chan- 
nel H . In contrast to the "confining phase" ran is not 
necessarily the mass of a meson. It corresponds to the 
eigenvalue of the ground state energy in the channel H. 
The conjecture made in ref. [10[ which was confirmed 
there can be translated to the following conjecture: When 
T/T c > 1 , in the "conformal region" defined by 



ma < cA 



IR 



(2) 



where c is a constant of order 1, the propagator Gn(t) 
behaves at large t as 



G H {t) 



_ exp(— rfiHt) 
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(3) 



which is a power-law corrected Yukawa-type decaying 
form instead of the exponential decaying form observed 
when rns > cAir,. 



culations in the N 



1 



2 case on a 16 x 64 lattice for 



several sets of parameters. The algorithm we employ is 
the blocked HMC algorithm [Hj]. 

After confirming the chiral transition occurs around 
K = 0.151 at /3 = 6.0, we choose the following values of 
P's: (3 = 6.5 (T - 2T C ), = 7.0 (T ~ 4T C ), (3 = 8.0 
(T ~ 16 T c ), (3 = 10.0 (T - 100 T c ), and /3 = 15.0 (T ~ 
10 5 T c ). We take several values of K for each /3 in such a 
way that the quark masses m's take values 0.00 < m q < 
0.10. The crude estimates of the temperature are based 
on A/? ~ 0.5 for the scale change of a factor 2 in the 
lowest perturbation. 

We choose run-parameters in such a way that the ac- 
ceptance of the global metropolis test is about 70%. The 
statistics arc 1,000 MD trajectories for thermalization 
and 1,000 MD trajectories or 500 MD trajectories for the 
measurement. We estimate the errors by the jack-knife 
method with a bin size corresponding to 100 HMC tra- 
jectories. 

We define the effective mass m//(i) by 



co8h(m H (t)(t - N t /2)) 



G H {t) 



coBh(m H (t)(t + l-N t /2)) G H {t+l) 

FIG. [2] shows the t dependence of the effective mass 
for the PS channel in the case T ~ 100T C with three 
types of sources. On the left panel, we see the clear 
plateau of the effective mass at t = 22 ~ 31 when quark 
mass is relatively large; m q = 0.30 (K = 0.125). On 
the other hand, on the right panel, we see the effective 
mass is slowly decreasing without no plateau up to t = 
31 when the quark mass is small; m q = 0.028 (K = 
0.135), suggesting the power- law correction. The power- 
law corrected fit for the local-local data in the K = 0.135 
case with the fitting range t = [15 : 31] with an = 1-1(1) 
reproduces the data very well. 

We show the effective mass plot for the local source - 
local sink in two more cases for T ~ 4T C and T ~ 10 5 T c 
in Fig- (J3J), together with the power-law corrected fit with 
the fitting range t = [15 : 31]. Both of the fit for f3 = 7.0 
at K = 0.142 with a H = 0.7(1) and the fit for (3 = 15.0 
at K = 0.130 with an = 1.21(4) reproduce the data for 
t = [15 : 31] very well. 

For all cases with 0.0 < m q < 0.10, we observe the 
similar power corrected behavior as shown in Fig. 2 
and we confirm the fit similar to those shown in Fig. 



Beta=10.0, K=0.125, Nf=2, 16 J x64, PS-channel 



Beta=10.0, K=0.135, Nf=2, 16 J x64, PS 
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FIG. 2. The effective mass for JV> = 2 at /3 = 10.0 (left: for if = 0.125; exponential decay) and (right: for K = 0.135; 
power-law corrected Yukawa decay): With three types of sources; the local-sink local-source (black squares), local-sink doubly- 
smeared-source (red circles) and local-sink doubly-wall-source (green triangles) 
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FIG. 3. The effective mass plots for local-sink local-source case and fits by power-law corrected Yukawa type decay (left: for 
f3 = 7.0 at K = 0.142) and (right: for /3 = 15 at K = 0.130) 



3 well reproduces the data. Thus we have verified our 
conjecture that the propagator on a finite lattice in the 
conformal region given by cq.(2) behaves at large t as 
a power-law corrected Yukawa-type decaying form given 
by cq.(3). We have not systematically estimated the con- 
stant c in eq.(2). In general, the c depends on Nf and 
P. We quote c ~ 2.0 x 2tt at /3 = 6.0 for Nf = 7 and 
(3=11.5 for N f = 16 [l(j for later use. 

Our numerical simulation confirms that the behavior 
of Gjy(i) is similar to each other between the small Nf 
QCD at high temperature and the large Nf QCD with 
an IR cutoff. One reason behind is that once gauge cou- 
pling constant is at the IR fixed point, the renormaliza- 
tion group equation is characterized by only one number, 
the mass anomalous dimension which is related to an- 
All the results show the existence of the " conformal re- 
gion" for a wide region for T > T c on a finite lattice. We 
believe it will survive in the thermodynamic limit from 
our theoretical argument. We hope to revisit more rigor- 
ous treatment of the thermodynamic limit in the future 
work. 



to point out that there is one difference between the small 
Nf QCD at high temperature and the large Nf QCD with 
an IR cutoff in the lattice size dependence because the 
origin of the IR cutoff is quite different. To understand 
this point, let us consider the phase diagram in terms of 
(3 and K for the small N f (N f = 2-6) QCD on a finite 
lattice. The chiral phase transition point on the massless 
quark line moves towards larger /3 with larger lattice, 
since the beta function is non-positive all through fi at 
T = T c . (Fo r Nf = 2 the chiral transition point j3 c t = 
3.9 - 4.0 [13 at N t = 4 moves to /3 ct ~ 6.0 at N t = 64. 
For N f = 6 j3 ct ~ 0.3(1) [13 at N t = 4 and the movement 
of (3 c t is very slow, which might be due to the smallness of 
the beta function as is in the walking technicolor theory.) 
Assume we are in the deconfining phase at some f3 with 
some lattice size. When the lattice size is increased with 
j3 fixed, the deconfining phase changes to the confining 
phase at some lattice size. This is clearly different from 
the large Nf QCD within the conformal window in which 
the IR fixed point remains even in the infinite volume 
limit at zero temperature. 



In relation to the thermodynamic limit, wc would like 



Now wc would like to discuss the implication of the 



existence of the IR fixed point in high temperature QCD 
and the existence of the "conformal region" on a finite 
lattice. In this article we concentrate on the issue of the 
order of the chiral phase transition in the Nf = 2 case. 
Besides our numerical results, our key observation is the 
existence of an IR fixed point at T > T c . We stress that 
the reasoning for the existence can be justified even in 
the thermodynamic limit. 

Pisarski and Wilczek [13[ mapped Nf = 2 QCD at high 
temperature to the three dimensional sigma model and 
pointed out that if £7^(1) symmetry is not recovered at 
the chiral transition temperature, the chiral phase tran- 
sition of QCD in the Nf = 2 case is 2nd order with 
exponents of the three dimensional 0(4) sigma model. 

The 0(4) scaling relation was first tested for staggered 
quarks [14j . For the Wilson quarks it was shown that the 
chiral condensate satisfies remarkably the 0(4) scaling 
relation, with the RG improved gauge action and the 
Wilson quark action j!5| and with the same gauge action 
and the clover-improved Wilson quark action [16| (See, 
for example, Fig. 6 in ref. [Iff). It was also shown for 
staggered quarks the scaling relation is satisfied in the 
Nf = 2 + 1 case [17] . extending the region from T/T c > 1 
adopted in [lf| and [Rf to the region including T/T c < 1 . 
These results imply the transition is second order. 

However, recently, it was shown that the expectation 
value of the chiral susceptibility x-n ~ Xa is zero [181 ] in 
thermodynamic limit when the SU(2) chiral symmetry 
is recovered under the assumptions we will discuss be- 
low. This is consistent with that the Ua{X) symmetry 
is recovered, which implies the transition is 1st order ac- 
cording to [l3| . Apparently the two conclusions are in 
contradiction. 

Here we revisit this issue with the new insight of con- 
formal field theories with an IR cutoff. We have to check 
two points: (1) whether the mass region investigated be- 
longs to the conformal region; (2) whether the power- 
corrected behavior of the propagators is properly taken 
into account. The chiral transition occurs at the mass- 
less quark line and therefore the mass region, where the 
0(4) scaling holds, will be in the conformal region. We 
took the quark mass region m q < 0.4 in the fit of Fig. 6 of 
ref. [lfl . On the other hand, the critical pion mass ra v for 
the boundary of the conformal region is estimated using 
the general equation m^ = c Air with Air = (N 3 xN t ) 1 / 4 
as in [Hj since the spatial lattice size is small 8 3 with 
N t = 4. It gives TOtt ~ 1.9, with c ~ 2.0 x 2n. Since the 
quark mass region m q < 0.4 corresponded to m^ < 1.7, 
we conclude that the mass region studied there belonged 
to the conformal region. This explains the fact that more 
than 30 data were excellently on the given smooth scaling 
curve with only one fitting parameter, the chiral transi- 
tion coupling constant, besides the scales of x-axis and 
y-axis. 

One more point is concerned with the power corrected 
Yukawa-type decay. We calculated the chiral condensate 



subtracted using 



ib = 2m q aZ 2 ; 

x,t 



[ir(x,t)ir(0)), 



where the sum is over space-time. This is equivalent to 
the time integration of the propagator with power cor- 
rected Yukawa-type decay. Note that this formula itself 
is always valid irrespective of the form of the propagator, 
but the scaling of the summed up quantity is affected by 
the power corrections. In [15j . we have directly studied 
the chiral condensate without noticing the power cor- 
rected decaying form, but the sum automatically takes 
into account the effects. Thus the analysis there is con- 
sistent with the new insight of conformal field theories 
with an IR cutoff and therefore the conclusion does not 
change that the 0(4) scaling relation is remarkably sat- 
isfied and the phase transition must be 2nd order. 

On the other hand, it is assumed in ref. [18| that the 
vacuum expectation value of mass-independent observ- 
able is an analytic function of mi, if the chiral symmetry 
is restored. However, in the conformal region the propa- 
gator of a meson behaves as eq.(3) and the relation be- 
tween the ran and the m q is given by the hyperscaling 
relation fl9|] 



m H 



,1/(1+7*) 



with 7* the anomalous mass dimension. This anoma- 
lous scaling implies m# is not analytic in terms of m 2 
and the analyticity assumption does not hold. It should 
be noted that the Ward -Takahashi Identities in [l8| are 
proved in the thermodynamic limit and therefore the nu- 
merical verification of the hyperscaling in the limit will 
be decisive. We stress however that the hyperscaling is 
theoretically derived with the condition of the existence 
of the IR fixed point and multiplicative renormalization 
of m q . We believe that this violation of the analytic- 
ity assumption resolves the apparent discrepancy as also 
mentioned in [18[ as a viable possibility. 

In the article in preparation we hope to discuss other 
physical implications of conformal field theories with an 
IR cutoff realized at high temperature QCD which are 
not discussed in this article, together with the discussion 
on the relation between large Nf QCD within the con- 
formal window and small Nf QCD at high temperature. 
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